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ABSTRACT 



We summarize the theory of coronagraphic optics, and identify a dimensionless fine- 
tuning parameter, which we use to describe the Lyot stop size in the natural units 
of the coronagraphic optical train and the observing wavelength. We then present 
simulations of coronagraphs matched to adaptive optics {AO) systems on the Calypso 
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1.2m, Palomar Hale 5m and Gemini 8m telescopes under various atmospheric condi- 
tions, and identify useful parameter ranges for AO coronagraphy on these telescopes. 
Our simulations employ a tapered, high-pass filter in spatial frequency space to mimic 
the action of adaptive wavefront correction. We test the validity of this representa- 
tion of AO correction by comparing our simulations with recent ii'-band data from 
the 241-channel Palomar Hale AO system and its dedicated PHARO science camera 
in coronagraphic mode. Our choice of monochromatic modeling enables us to distin- 
guish between underlying halo suppression and bright Airy ring suppression in the final 
coronagraphic images. For a given telescope-^0 system combination, we find that AO 
systems delivering images with Strehl ratios below a threshold value are not well-suited 
to diffraction-limited coronagraphs. When Strehl ratios are above this threshold, an op- 
timized coronagraph with occulting image plane stops as small as A\/D create a region 
around the AO target where dynamic range is significantly enhanced. 

Subject headings: instruments: coronagraphs - instruments: miscellaneous - tech- 
niques: coronagraphy 

1. Introduction 

Many ground-based telescopes use adaptive optics (AO) to produce diffraction-limited images 
at near-infrared and visible wavelengths. The point spread function (PSF) of a telescope using 
AO typically consists of a bright, diffraction- limited core with several Airy rings superimposed on 

a wide scattered light halo containing several percent of the total flux. The improvement in image 
quality that AO provides makes it possible to study the region within a few times the diffraction 
width of the image of a bright star, with dynamic range limited by the presence of the halo and 
bright Airy rings rather than by atmospheric seeing. 

Present-day astronomical AO systems routinely deliver Strehl ratios of 30-70% (the Strehl ratio 
is the ratio of the peak intensity of the image to the peak intensity of the image if the wavefront 
were free of all aberration) . In this paper we look at what is possible when optimized coronagraphs 
are used with AO systems delivering 50-95% Strehl ratios. We believe that this should be the 
next scientifically-driven step for traditional AO systems, since the instruments we model here will 
open up a new range of ground-based astronomical investigation, and provide a stepping stone to 
even higher dynamic range ground-based astronomy that is outside the reach of multi-conjugate 
AO systems used to correct wide fields of view (Berkefeld et al. 2001a). 

A coronagraph used in conjunction with an AO system can improve the sensitivity of an 
imaging system to faint structure surrounding a bright source. This device blocks the core of the 
image of an on-axis point source and suppresses the bright diffraction rings and halo, removing 
light which would otherwise reduce the dynamic range of the imaging. This enables faint off-axis 
structure to be observed. 
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Initial observational investigation of the dynamic range achievable with a coronagraph on the 
Palomar Hale 5m AO system has demonstrated the importance of understanding the interaction 
between the AO system and coronagraphic instrument parameters (Oppenheimer ct al. 2000). 
The choice of pupil plane and image plane stops sizes can only be made by understanding image 
formation by the AO system and the interaction between the AO system and the coronagraph. 
A traditional coronagraph reduces off-axis light from an on-axis source with two optical stops: an 
occulting stop in the image plane, and a matched Lyot stop in the next pupil plane citeplyot. An 
alternative solution to this problem is a nulling coronagraph (e.g., Guyon et al. (1999)), although 
this novel approach still needs technological development in order to be applied to broad-band 
imaging. The discovery space of a nulling coronagraph is also somewhat complementary to the 
optimized AO coronagraphs we describe here. 

With an AO system in place, the image plane stop can be made very small: only slightly larger 
than the diffraction spot itself. However, a pupil plane stop matched to a reduced image plane stop 
must have a small clear diameter to significantly reduce the off-axis throughput; there is a trade- 
off between throughput and scattered light suppression. In some situations using a coronagraph 
does not improve the final dynamic range of an instrument — it simply extends the exposure time 
required for a given amount of detected signal. 

To address this issue, we simulated image formation in a coronagraph mounted on a telescope 
with an AO system, and investigated the effects of atmospheric turbulence under a range of seeing 
conditions, telescope sizes, and AO system performance levels. We introduce the principles of 
coronagraphy with a simple one-dimensional model, and discuss the results of the simulations in 
the context of current and next-generation instruments. 

Our model is an extension of that of Nakajima (1994), who modeled low order ground-based 
AO coronagraphs to determine the detectability of Jupiter-mass companions around main sequence 
stars. Nakajima simulated the AO wavefront correction by zeroing out the lowest order Zernike 

coefficients in the expression used to generate realizations of Kolmogorov-spectrum fluctuations in 
the atmosphere's refractive index. In Nakajima's work, the Lyot stop is "oversized" by a fixed 
amount (10%), and occulting stop sizes are chosen to be 1, 5, 10 and 15 times the telescope 
resolution. We extend Nakajima's approach in four ways. First, we match the Lyot stop oversizing 
to the size of the image plane stop to optimize coronagraphic performance. Second, we use a graded 
high-pass filter to mimic AO correction. Third, we investigate higher order adaptive correction of 
the incoming wavefront to model AO systems with a few thousand actuators. We concentrate on 
the smallest stop sizes allowed by Fourier optics that do not lead to an unacceptable reduction in 
throughput. In earlier work (Makidon et al. 2000) we investigated occulting stop sizes between 
3X/D and 6X/D radians. Here we present results for a single occulting image plane stop with a 
diameter of AX/D radians. Present-day computer memory and processor speeds, as well as the 
development of a rapid Markov algorithm to create the Kolmogorov-spectrum phase screens used 
to simulate atmospheric effects (Glindcmann ct al. 1993; Berkcfcld et al. 2001b) enabled us to 
generate 1000 realizations of each instrument configuration. We therefore simulate long exposure 
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AO images which can be used to estimate detection hmits, since speckle noise is often the hmiting 
factor in faint companion detection (Racine et al. 1999; Marois et al. 2000). 

2. Theory 

To set the stage for discussion of a more reahstic instrument, we illustrate the Fourier optics 
of a one-dimensional coronagraph. Other more formal expositions of coronagraphic imaging can be 
found in the literature (Noll 1976; Wang & Vaughan 1988; Malbet 1996). Our analysis assumes 
that the Praunhofer approximation applies, i.e., that the transverse electric field in the image plane 
is the Fourier transform of the phasor of the wavefront phase in the pupil plane (if (j){x,y) is the 
phase, then e"^^^'^^ is the corresponding phasor). We make use of the standard Fourier analysis 
results which can be found in Bracewell (1986). 

For this analysis we consider only monochromatic imaging, but note that in typical broad-band 
imaging, the final image can be described by the sum or integral of several monochromatic images, 
weighted by the instrumental transmission function. Wavelength variation across the band will act 
in such a way as to smear image features radially by the same factor as the fractional bandwidth, 
since the wavelength enters into diffraction-limited image formation only in the combination (X/D). 
As a result, bright Airy rings will get wider, but coronagraphic suppression of such rings will persist. 
By treating the monochromatic case we can distinguish clearly between halo suppression and Airy 
ring suppression. Secondary support spiders and scintillation (field strength variation) are not 
modeled here. 



2.1. The one-dimensional coronagraph without an atmosphere 

In the absence of atmospheric degradation, a monochromatic on-axis source at infinity produces 
a transverse electric field at the telescope pupil 

E = Eoi?e(e^('=^-'^*)), (2-1) 

where k = 2tt/\, the z axis is the optical axis, u the angular frequency of the wave, t the time, and 
ReQ denotes the real part of a complex number. 

We follow the passage of the incident wave's field through a one-dimensional coronagraph. 
Figure 1 shows a diagram of the optical path. 

We label eight key locations along this path with letters a-h. Eight plots in the figure show 
the electric field due to an on-axis source at these key locations, and the transmission functions of 
the optical stops that affect the incident wave as it passes through the coronagraph. 

First, the incoming wave passes through the telescope aperture (Figure la). We represent this 



-5- 



interaction by multiplying the field by the aperture stop function, so that in the pupil plane 

„ ,^ , , n(x) = 1 for Ixl < 1/2, 

= ^n(x/Z>.), where ^J^j ^ ^ (2-2) 

Here Dx = D/\ is the number of wavelengths across the telescope aperture (which is also the 
inverse of the angular resolution of the telescope as measured in radians) . We denote a pupil plane 
coordinate by x, and an image plane coordinate by 0. The telescope optics then form the wave into 
an image (Figure lb). The electric field in the image plane is the Fourier transform of the aperture 
field Eai 

Ef, oc smc{DxO), (2-3) 

where 9 is the field angle in radians in the first image plane. We omit the constants of proportionality 
for simplicity. 

In a conventional imaging camera this image field would fall on a detector here. However, in 
a coronagraph the star is occulted by a field stop in this image plane. We describe the stop in 
terms of a shape function w{Dx6/s), which is unity where the stop is opaque and zero where the 
stop is absent. If w{6) has a width of order unity, the stop size will be of the order of s resolution 
elements. The transmission function in the image plane is therefore 1 — w{Dx9/s) (Figure Ic). To 
illustrate the present discussion, we take w{6) = exp(— ^^/2). The field in the first imaging plane 
after the occulting stop (Figure Id) can therefore be written as 



Ea oc smc{Dxe){l - w{Dxe/s)). (2-4) 

This occulted image is relayed to a detector through a second pupil plane. The electric field 
at this second pupil is the Fourier transform of the occulted image field (see Figure le): 



E, oc Il{x/Dx) * {5{x) - W{s x/Dx)) (2-5) 

Here W is the Fourier transform of the image stop function 6{x) is the Dirac delta function, 
and * denotes convolution. W has width of order unity, although it will not have bounded support 
for occulting stop shape functions of finite extent, such as the hard-edged stop w^O) = H(0) (the 
support of a function is the set of points at which the function is non-zero). The geometrical 
significance of (2-5) becomes clear if we rewrite it as 



Ee OC U{x/Dx) - Il{x/Dx) * W{s x). (2-6) 

If the image stop is completely opaque at its center, w{0) = 1. This means that its transform, 
(s/Dx) W{sx/Dx), has unit area, regardless of any re-scaling of the argument of w. This makes 
for cancellation of the field across most of the pupil when s » 1. In Figure 2 we show how the 



-6- 



equation 2-6 is constructed graphically, using a Gaussian image stop whose width is bX/D {i.e. 
s = 5). This shows why the Lyot stop must mask out a border of order D/ s wide around the pupil 
boundary to produce significant reduction in the throughput of unoccultcd light from the on-axis 
source. It is only at this stage that the coronagraph increases the dynamic range of the final image. 

The unoccultcd light (Figure Id) has a highly periodic distribution, with periodicity ~ \/D. 
In the following pupil plane (which is the transform space of the image plane), this energy is 
concentrated near ibD/2. The larger the occulting stop diameter, the more E^i looks like a pure 
sinusoid, and the more the unocculted energy is localized in the neighborhood of the boundary of 
the following pupil. 

In seeing-limited coronagraphs, the occulting stop is typically many diffraction widths in size 
(s > 10). Consequently, the Lyot stop need only be undersized by a small fraction of the pupil 
diameter {e.g. 10% or less), resulting in minimal loss of throughput for unocculted, off-axis sources. 
In contrast, off-axis throughput in an optimized, diffraction-limited coronagraph with significant 
rejection of on-axis light must fall dramatically as the image plane stop shrinks to a few diffraction 
widths. This is because the spillover of unocculted on-axis light occurs in a wide border around the 
pupil boundary in the plane of the Lyot stop, so the Lyot stop must obscure a sizeable fraction of 
the re-imaged primary mirror to remove the on-axis spillover. Hence, off-axis throughput is reduced 
as well. 

Since the scale of the Lyot stop oversizing is D/s, we fine-tune the Lyot stop diameter so 
that it obscures a border TD/s around the perimeter of the primary. The Lyot stop diameter is 
therefore 

DLyof = D- 2TDls. (2-7) 

If a secondary obstruction is present, then the Lyot stop must block out a similar border around 
the inside edge of the annular pupil. This is why small secondary mirrors benefit diffraction-limited 
coronagraphy. In section 4.1 we describe how to choose an optimum value of T for a given telescope 
and AO system under given atmospheric conditions. 

Using the theory outlined above, in the the case of an unobstructed primary aperture with an 
image plane stop of 5A/D {i.e., s = 5), approximately 16/25 of the aperture should be obscured 
by a matched Lyot stop: when projected back onto the primary pupil, the Lyot stop is opaque 
outside a circle of diameter ~ 3D/5. Rejection of unwanted on-axis light must be balanced by 
signal-to-noise considerations pertaining to the off-axis source brightness. This places a practical 
lower limit on the angular size of the occulting spot in the first image plane. 

The above arguments hold for two-dimensional apertures as well. The derivation is analogous 
to the one-dimensional case, although the functions and transforms become two-dimensional {e.g., 
for a circular telescope aperture, the sine function is replaced by the Airy function). Wang & 
Vaughan (1988) describe the two-dimensional case, and Malbet (1996) treats the PSF's of off- 
axis sources in such coronagraphs. 



- 7- 



2.2. The coronagraph in the presence of an atmosphere 

In the absence of scintillation, we can model atmospheric effects on an incident plane wave by 
multiplying the aperture illumination function n(x) by an atmospheric disturbance phasor e '"^^^^ . In 
our one-dimensional example, the resultant image field is the Fourier transform of e^'^^^^Il{x/Dx)- 
For telescope diameters larger than ~ 0.5m at wavelengths shorter than ~ 2 microns, the atmo- 
spheric phase function typically has significant power at length scales smaller than the aperture 
size, so a long exposure image under these circumstances exhibits the familiar unimodal seeing 
disk that is typical of images from large ground-based telescopes. In consequence, coronagraphs 
without AO systems on seeing-limited telescopes use occulting stops many diffraction widths wide, 
with concomitant high throughput as explained earlier. Stellar coronagraphy without adaptive 
wavefront correction, or with only tip-tilt systems, is still useful to prevent detector saturation, and 
has produced observations of considerable scientific value (Golimowski et al. 1992; Nakajima et 
al. 1995). 

Ideally, a coronagraph with an AO system would be just like the coronagraph without an 
atmosphere. The image of a point source in the first image plane would be a pure Airy disk, and 
the size of the Lyot stop could be chosen using simple Fourier theory and a diagram like Figure 
1. However, an AO system cannot correct atmospheric effects on spatial scales smaller than the 
inter-actuator spacing in the pupil plane. Power at high spatial frequencies that goes uncorrected 
in the pupil plane transforms into noise on large angular scales in the image plane. The corrected 
image consists of a diffraction limited core surrounded by an extended halo due to the uncorrected 
aberrations, a remnant of the uncorrected seeing-disk. The size and shape of the halo reflect the 
number of actuators in the AO system and the moment-to-moment characteristics of the turbulence 
in the atmosphere. 

In the presence of this halo the choice of Lyot stop is not obvious. Indeed, a given AO 
configuration may not be able to reduce the power in the extended halo sufficiently to justify 
diffraction-limited coronagraphy — coronagraphy where the occulting mask is only a few diffraction 
widths across. In order to predict the size of the uncorrected halo and to understand how it affects 
coronagraphy, we modeled the effects of atmospheric turbulence on the incoming wavefront and the 
correction of these aberrations by AO. Our models allow us to find useful operating parameters for 
diffraction-limited coronagraphs working in concert with AO systems for a range of telescope sizes, 
AO capabilities, and observing conditions. Sections 4.1 and 4.2 discuss these issues in some detail. 

3. Numerical Simulations 

Each of our numerical simulations is characterized by the telescope entrance pupil size and 
geometry, the seeing, D/vo, the number of actuators across the primary, Nact, the size of the 
occulting image plane stop, s, the Lyot stop fine tuning factor, J^, and the linear size of the array, 
Ng, sampling the incoming wavefront. We generate 1000 independent realizations of Kolmogorov- 
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spectrum phase screens in 2Ns by 2Ns arrays, with the spatial samphng of these arrays chosen to 
provide several samples across each Fried length Vq (Fried 1966), while constraining the sampling 
to force Ng to be a power of 2. Our resultant spatial sampling ranges from 3.2 to 21 samples across 
To (see Table 1). 

We Fourier transform the input phase arrays, and multiply them by a high-pass filter to 
mimic the action of AO. Then we reverse-transform the filtered arrays to obtain the AO-corrected 
wavefront. 

The shape of the AO filter near cutoff depends on details of the AO system. Deformable 
mirror (DM) actuator infiuence functions which extend to neighboring actuators' positions reduce 
the sharpness of the cutoff. Noisy wavefront sensing and intrinsic photon noise reduce the efficacy 
of high spatial frequency wavefront correction. The flow of the atmosphere past the telescope pupil 
and a non-isotropic refractive index spatio-temporal distribution also change the shape of the AO 
filter, as does imperfect DM calibration. 

The cutoff frequency of the AO filter cannot be higher than the spatial Nyquist frequency of 
the actuator spacing, kAO = ^act/^D. At a given wavelength A, this spacing corresponds to an 
angle 

Oao = NactX/2D. (3-1) 

In earlier work (Makidon et al. 2000), our AO high-pass filter was the complement of a Hanning 
filter. This has a continuous derivative everywhere, and a very smooth approach to the cutoff. 
As we show in Section 5, this proved to be too conservative, as it underestimates the amount of 
observed AO correction. In these simulations we use a parabolic filter: 

A{k) = {k/kAo)"^ for k < kAO 
A{k) = 1 otherwise. 

We show that this matches observed results over the range where comparison with data is valid. 

To avoid edge effects introduced by the Fourier filtering, only the central Ng by A^^ section of 
the filtered array is retained. This is then multiplied by a binary mask representing the telescope 
entrance pupil. This mask is opaque (zero) outside the primary mirror edge, and inside the sec- 
ondary obstruction. Secondary support spiders, mirror surface micro-roughness and scattering are 
not considered in these simulations. 

We embed the filtered and masked array in the center of an SA^ by BA^ zero-filled array. This 
results in image field sampling of X/8D. which allows for effective comparison with data taken with 
^ X/4D pixels. 

We then create the complex phasor e"^^^'*') describing the electric field corresponding to the 
phase <p{x, y), and Fourier transform the phasor to obtain the first image field. We average all 1000 
realizations of the intensity of this image, as well as the Strehl ratios derived from these images as 
a measure the AO system performance. 



(3-2) 
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This ^O-corrected image field is then multiplied by another binary mask to produce the effect 
of the image plane stop, and the product is inverse-transformed to create the pupil field at the 
Lyot plane, where it is multiplied by a third binary mask representing the Lyot stop. The final 
image is produced by Fourier transforming this field. We average all realizations of the final image 
intensity, along with the values of various throughput and geometrical obstruction descriptors of 
the optical system. These averaged images can be thought of as a single image with an exposure 
time corresponding to a thousand times the speckle lifetime. This is equivalent to a ~ 100 second 
exposure at Palomar in the ii'-band. 

4. Optimizing Coronagraph Design 

In this section we present the results of numerical investigations into when AO benefits coron- 
agraphy, and how one optimizes the design parameters for a coronagraph used in conjunction with 
an AO system. 

4.1. Optimizing the Lyot Stop 

In section 2 we defined the diameter of the Lyot stop as D — 2TD/s, where T is the Lyot stop 
fine tuning factor. Here we show how the value of is chosen to optimize a coronagraph for a 
specific telescope under particular atmospheric conditions. We present the results of an exploration 
of a coronagraph on the Gemini telescope with atmospheric turbulence described by D/ro = 30. 

The uppermost curve in Figure 3a shows an azimuthally averaged PSF that an AO system 
with 51 actuators across the primary diameter (a total of 2042 actuators) delivers. A hundred 

independent atmospheric phase disturbance realizations have been averaged into this PSF, which 
is normalized to unity at its center. The Strehl ratio of the corresponding image is 82%. 

Below this we show azimuthally averaged PSF's of a coronagraph with a central stop of 4X/D 
diameter, with three different values of J^. These curves have been re- normalized to take into 
account the dimming of the image due to throughput losses that result from undersizing the Lyot 
stop. This is accomplished by dividing the simulated profile by the Lyot stop clear fraction (the 
ratio of the clear area of the Lyot stop, when projected back to the primary, to the clear area of 
the entrance pupil). If the AO-corrected PSF and a coronagraphic PSF coincide at some particular 
value of the radial distance from the central object, then, at that separation there is no dynamic 
range gain to be had from using that coronagraph. At such a location the same companion signal- 
to-noise ratio will be achieved when the same number of source photons have been detected. In 
this case the coronagraphic configuration simply increases the exposure time required to achieve 
the same signal-to-noise ratio, because it is equivalent to using a telescope with a smaller collecting 
area. 
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Figure 3b shows the ratios of the various renormahzed, azimuthally averaged coronagraphic 
PSF's to the PSF with just the AO system and no coronagraph. Three values of the fine tuning 
parameter are presented: 0.25,0.5 and 0.75. The PSF for least aggressive Lyot stop (J- = 0.25) 
shows the least dynamic range benefit. Bright re-imaged pupil edges are not sufficiently masked 
out, though the Lyot stop is 72% clear. The most aggressive stop = 0.75) masks out much more 
of the bright pupil edges, but the Lyot stop is only 20% clear. This results in a drastic reduction 
in the brightness of off-axis sources. The intermediate Lyot stop (.F = 0.5) does a better job of 
balancing the obscuring of bright pupil edges, and has a Lyot stop throughput of Lyot stop (43%). 
For a specific instrument project this optimization will need to be done on broad-band images, with 
finer resolution in the J- parameter space. We present this coarse optimization search as a model 
for further work tailored to particular telescopes and science drivers. 



4.2. Threshold AO performance for coronagraphy 

Here we consider the transition between seeing-limited and high Strehl imaging. We select fixed 
seeing conditions {D/ro = 30), and vary the actuator spacing of the AO system. The resulting 
images have Strehl ratios between 30% and 90% as the number of actuators across the primary 
diameter goes from 26 to 71 (Figure 4). The validity of our AO correction algorithm is discussed 
in section 5, we merely note here that it is based on real data. We normalize the PSFs in Figure 
4 to be unity at the center because we wish to present dynamic range gains rather than absolute 
flux levels. 

The PSF is described by a diffraction-limited core, a flat AO-corrected halo out to a shoulder 
at a radius 6ao where the AO stops working. After that the profile shows the familiar atmospheric 
wings. The fiat plateau of AO correction drops in intensity and grows in radius as the number of 
actuators increases. It is the uncorrected light in this area that is redistributed into Airy rings by 
the AO system, and it is this light that is not removed by a coronagraph. The image with a Strehl 
ratio of 30% has a ten-fold decrease in the intensity of this scattered light as compared with the 
53% Strehl image. This demonstrates why there is a threshold Strehl required for coronagraphic 
imaging targeted at a particular angular distance from the central source. 

Azimuthally averaged PSF's (dotted lines) for the same 100 realizations of the atmospheric 
phase disturbances corrected by AO systems with a range of actuator spacings are shown in Figures 
5a - 5f. Below these PSF's we plot renormahzed PSF's for a coronagraph with a AX/D occulting 
spot and a Lyot stop fine-tuning parameter of ^ = 0.5. The vertical scales on all these plots is the 
same, and the PSF is normalized to unity at the origin. 

The effect of increased AO correction on direct imaging is manifest in the steadily decreasing 
intensity of the wings relative to the central peak. 

At a Strehl ratio of 53% (Figure 5a), the wings of both the direct and coronagraphic PSF's 
are about 10"^ of the central intensity. There is little suppression of the halo of the image, though 
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bright Airy rings are removed out to a radial distance of ^ 6X/D. Scattered light from uncorrected 
atmospheric aberrations are not suppressed by the coronagraph. At a Strehl ratio of 62% (Figure 
5b), AO correction produces visible Airy rings further out, since more light is guided into the Airy 
pattern and less light is scattered semi-randomly into the halo. Again, residual scattered light is 
not suppressed coronagraphically, and it still limits dynamic range. In Figures 5c and 5d we see 
further reduction in the scattered light halo, and a corresponding increase of energy in the pure 
Airy pattern. Coronagraphic imaging suppresses the image to well below the darkest Airy rings, 
a phenomenon we call halo suppression, since that is how it appears to the observer, even though 
it is really suppression of a coherent diffraction pattern emerging because of the reduction in halo 
intensity. 

By the time the Strehl ratio has risen to 82% (Figure 5d), an annulus of greatly increased 
dynamic range opens up between the occulting stop edge at a radius of 2X/D, and the radius at 
which the remaining halo is comparable in intensity to the Airy pattern, at a radius of AX/D. It is in 
such regimes that dynamic range increases dramatically. Below such Strehl ratios a coronagraphic 
spot with a AX/D diameter would not produce much gain in dynamic range. 

The improved dynamic range would continue indefinitely if AO guide stars were bright enough. 
However, for a given primary area, as the number of AO channels increases, the number of photons 
per channel available to use for AO correction decreases. We extrapolate from Palomar AO system 
performance to the regimes shown in Figures 5e and 5f as goals that coTild be reached by the next 
generation of AO systems on 8m telescopes at good astronomical sites. 

5. Comparison of Simulations with Palomar PHARO Observations 

We compared our monochromatic models to a 180 second i^-band exposure obtained with 
the PHARO camera on the Palomar AO system, using an image plane stop of (X'96 in diameter 
(Oppenheimer et al. 2000). We used the appropriate Lyot oversizing to match the PHARO coro- 
nagraphic Lyot stop dimensions {!F = 1.07, or a clear aperture of 4m), and match the AO system 
actuator spacing of 16 across the primary. Details of the PHARO camera optics are described in 
Hayward et al. (2000). 

Figure 6 shows an overlay of the radial profiles of our simulations with D/rg = 10 and the 
parabolic AO high pass filter (dashed), against the data (solid). Both images were normalized to 
contain unity power out to a radius of 4'.'0. 

Our choice of D results in a match with the wings of the stellar image past 9ao = 8X/D 
(where the AO ceases to improve the image). We did not need to fine-tune the Fried length 
to match the observing conditions when the data were taken, as our initial guess at provided a 
sufficiently accurate estimate of the azimuthally averaged stellar wing profile. Other effects (non- 
Kolmogorov spectra, instrumental scatter, the effects of the secondary support spiders, waffle or 
other wavefront reconstructor errors, etc.) could also be present in the PHARO data. 
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The simulated image plane stop has a diameter of 10.6X/D at 2.2 microns. Therefore the 
region where our model is tested by the data lies between 5.3\/D ((y.'48) and 9ao = 8X/D ((/.'TS). 
Our model differs from observations by a few percent over the entire range where the compari- 
son is valid. Detailed work on matching simulations with this data will be presented elsewhere 
(Sivaramakrishnan and Oppenheimer 2001). 

The comparison with data suggests that our representation of AO correction as a parabolic 
filter is accurate enough to be used to predict system performance of the next generation of AO 
coronagraphs. 

6. DISCUSSION 

In order to span a range of telescope sizes, we simulated the KPNO Calypso 1.2m, the Palomar 
Hale 5m, and the Gemini 8.1m. We reproduce the secondary obscuration ratios of each telescope, 
though we do not model secondary support spiders. 

Atmospheric phase screen realizations using one or a few representative values of the Pried 
length relevant to the sites are input to the optical train simulation after being smoothed by the 
AO filter. The effect of wind on the delivered Strehl ratio have not been incorporated. 

We match the AO system actuator density on the primary apertures to those presented in 
previous work (Makidon et al. 2000). A summary of the relevant input parameters (primary and 
secondary mirror diameters, number of actuators, and input D/vq values) can be found in Table 1. 

While we examined a variety of image plane occulting stop sizes, wc onh^ present one stop 
size, AX/ D, in this paper. Our choice of this stop diameter is based on previous work which 
suggests AX/D is close to the smallest effective stop size for imaging stellar coronagraphs with 
Strehl ratios between 70% and 95% (Makidon et al. 2000). Larger stops simply obscure more of 
the improved image without a significant increase in the dynamic range of the final images. We find 
that stops smaller than AX/D require extreme oversizing of the Lyot stop, dramatically reducing 
system throughput. We note that for larger values of the image plane stop diameter, the Lyot 
throughput penalty is reduced. In consequence, the optimal value of T is expected to grow with 
increasing values of s when s < 10. 

Two values of 0.5 and 0.25, are examined here. These yield Lyot stop outer diameters of 
0.75 D and 0.875 D respectively for our image plane stop. The effect of the two Lyot stops on the 
final image are illustrated in the PSF plots of Pigures 7, 8 and 9. In general, we find a coronagraph 
coupled to the more aggressive Lyot stop {J^ = 0.5) is significantly better at suppressing both 
Airy ring and halo contributions to azimuthally averaged PSFs than the higher throughput one 
(:F = 0.25). 

The coronagraphic PSF^s presented in Pigures 7, 8, and 9 take into account the Lyot stop 
throughput losses, as they have been renormalized by dividing by the ratio of the Lyot stop area to 
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the primary aperture area. When the image profiles with and without the coronagraph coincide, 
no gain in dynamic range is attained by using a coronagraph. 

We indicate the stop edge in these figures with a vertical line. We also mark where the AO 
stops working (at the angular distance 9ao) with a triangular symbol if this occurs within the plot 
boundaries. In our Gemini simulations AO this angle falls outside the borders of the PSF plots. 

6.1. Kitt Peak Calypso 1.2m telescope 

The Kitt Peak Calypso 1.2m telescope was designed and constructed to minimize the dele- 
terious effects of ground-level turbulence and scattered light on astronomical observations (Smith 
1999). This telescope's instrumentation operates in optical bandpasses. 

For a small-aperture telescope like Calypso, we found that Strehl ratios of at least 60% are 
necessary for useful AO coronagraphy. The simulation we present has a 75% Strehl ratio. Smith 
(1999) reports upper quartile ro's of 20cm or less in the V band at this site — we use a value of 
24cm. 

The Airy rings are well-suppressed out to many diffraction widths when = 0.5 (Figure 7a). 
The broad underlying halo is reduced by about a factor of three a few resolution elements past the 
image plane stop with this Lyot stop. 

6.2. Palomar Hale 5m telescope 

The Hale 5m has a functioning 241-channel AO system with a dedicated near-infrared science 
camera, PHARO (Hayward et al. 2000), which contains a variety of coronagraphic image plane 
and Lyot stops. 

When D/rQ = 10 we expect Strehl ratios of the order of 80% (Figures 8a and 8b). Such 
conditions can be expected in the i^-band on good nights. At this Strehl ratio the Airy ring 
structure is almost completely removed, and halo is reduced by about a factor of two between the 
stop edge at 2X/D and ~ 5\/D when = 0.5 (Figure 8b). The benefits of using a coronagraph 
when the AO delivers a 50% Strehl ratio are greatly diminished, as can be seen by comparing figures 
8b and 8d. Thus the current 241-channel Palomar AO system is suited to -fC-band coronagraphy 
under good conditions. A system upgrade to e.g., twice the number of actuators will enable the 
coronagraph to open up a larger search annulus in the i^-band, and produce useful dynamic range 
gains in the if-band. To exploit the full potential of coronagraphic imaging on this system, upgraded 
wavefront reconstruction may be needed to delve into the cores of the images at angular separations 
of the order of the stop radius of 2A/-D from the target AO star. 
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6.3. Gemini 8.1m telescope 

Gemini has a near-infrared coronagraph under construction, and will be equipped with a low 
order AO system, which enables high order second stage AO (Angel 1992) to be incorporated 
into a dedicated AO coronagraph. We used values of the Pried length which apply in the near- 
infrared. Under these conditions, expected Strehl ratios of 80% - 95% produce extremely effective 
coronagraphic suppression. In addition, 9ao = 25X/D, so the annulus of high dynamic range 
imaging is relatively large (Figure 9a through 9f). 

With D/vo = 30, image halo reduction by a factor of three is observed a few diffraction widths 
outside the stop edge, and Airy ring suppression by another factor of three extends out to ^ 5X/D 
(Figure 9f). 

At the other extreme, when D/vq = 15 (e.g., in the X-band), performance improves dramati- 
cally — in Figure 9b we see Airy ring suppression of a factor of over 10 or more in the bright rings 
located at 3.5 and 5.5 X/D, and halo suppression of three out to ~ 7X/D. 

In Figures 10a - lOd we show simulated images for companions with brightness differences Am 
of -|-5, +7.5 and +10 magnitudes relative to the on-axis target. In the absence of noise, we see that 
a companion with Am = +10 at a separation of 0'.'16 is observable in the i?— band under good 
seeing conditions {D/vo = 15) with a 2000 channel AO coronagraph with a 0'.'16 diameter field 
stop (Figure lOd). The exposure time corresponding to these simulated images is ^ 1000 speckle 
lifetimes. 

7. Conclusion 

Diffraction-limited AO coronagraphy targets an entirely different search space from seeing- 
limited coronagraphs. Traditional coronagraphy under seeing-limited conditions provides image 
suppression where the seeing halo drops below the Airy pattern's wings (Macintosh 2000). Typ- 
ically this is tens of diffraction widths away from the central bright object. The advantage of 
diffraction-limited coronagraphic imaging is concentrated in an annulus starting just outside the 
image stop, and ending where the AO system stops improving the wings of the AO target star's 
image. The size of the outer edge of this annulus is set by observing wavelength and the spatial 
Nyquist frequency of the AO system actuators when projected on to the primary. When small 
image plane stops are used the Lyot stop size must be carefully chosen to find the best trade-off 
between throughput and image suppression. Coronagraphic instruments might work extremely well 
in the near-infrared and longer optical bandpasses on large telescopes with the next generation of 
deformable mirrors using a few thousand actuators, provided wavefront correction is done well. 

For 1-m class telescopes at IjUm, we find that AO systems can deliver Strehl ratios > 60% 
without excessively compromising the brightness limits of available AO . However, only modest 
contrast ratio gains are achieved with a coronagraph on such telescopes. Coronagraphic suppression 
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is evident mostly in the bright Airy rings of the image. 

The existing 241-actuator Palomar AO system is useful for iC-band coronagraphy. Stop sizes 
down to at least AX/D are theoretically possible. Improved wavefront reconstruction may be neces- 
sary for these small occulting stops. With double the number of AO channels and reduced wavefront 
sensor read noise, this telescope could deliver Strehl ratios of 90% at 2.2;um. H- and i^-band coro- 
nagraphy would then become more scientifically profitable because of the increased dynamic range 
in this search space. 

On 8-m class telescopes, tremendous coronagraphic gains are to be had with 2000 channel AO 
systems. With occulting stop diameters of ~ O'.'IS, three-fold halo suppression and ten-fold bright 
ring suppression can be expected a few diffraction widths from the stop edge. Such a system on 
the Gemini telescopes would be ideally suited to high dynamic range coronagraphy in the iJ-band, 
with limiting Strehl ratios of ~90% possible around an MO dwarf AO target ~20 parsecs away. 

We thank B. Brandl and S. R. Kulkarni for posing questions that led to this investigation, and 
P. E. Hodge, T. Nakajima, B. A. Macintosh, and B. R. Oppcnhcimcr for insightful discussions. B. 
R. Oppenheimer also provided us with pre-publication data. We are indebted to J. P. Lloyd for a 
critical reading of our simulation code, and the anonymous referee for helpful suggestions. We also 
thank R. J. Allen, J. R. Graham, E. P. Nelan, and R. L. White for encouraging the investigation. 
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Table 1. Properties of Telescope Simulations 



Tck-scopo 


D 


D,./D 


Nan 


-^^ rhd/i 




D/ro 


Sl.rohl (%) 


Calypso 


1.2 


0.36 


7 


38 


64 


5.0 


75 


Palomar 


5.08 


0.33 


16 


201 


64 


10.0 


79 














20.0 


48 


Gemini 


8.1 


0.13 


51 


2042 


128 


15.0 


94 














22.5 


89 














30.0 


80 



Note. — For this table, D and Dg are the primary and secondary mirror 
diameters (in meters). Nact is the number of actuators across primary mirror 
diameter, Nchan the estimated number of AO channels within the pupil, and 
Ng the number of wavefront phase samples across the primary mirror diameter. 
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Fig. 1. — One-dimensional coronagraph summary, with locations and electric field or stop profiles 
of: (a) primary pupil for on-axis source; (b) image before image plane stop; (c) image plane stop; 
(d) image after image plane stop; (e) pupil before Lyot stop; (f) Lyot stop; (g) pupil after Lyot 
stop; and (h) final on-axis image. In this example, 98% of the incident power is blocked by the 
coronagraph. 
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Fig. 2. — Graphical representation of the Lyot plane electric field calculation: (a) pupil function of 
width D; (b) Gaussian profile image plane stop with 5X/D standard deviation produces a Gaussian 
with standard deviation D/5 in the Lyot plane; (c) the convolution of the pupil function with the 
transform of the stop profile; (d) the final Lyot stop field showing bright edges and no energy in 
the center. 
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Fig. 3. — Optimizing the Lyot stop diameter. Top panel: azimuthally averaged AO-only PSF (dot- 
ted line), and three coronagraphic PSF's with different Lyot stop diameters. The coronagraphic 
PSF's have been re-normalized to take Lyot stop losses into account (see text). A A\/D angular 
diameter occulting image plane stop, whose edge is indicated by a dashed vertical line in both pan- 
els, is used in these coronagraphic simulations. Three values of the Lyot stop fine tuning parameter 
describing the Lyot stop size are presented. Bottom panel: ratios of the azimuthally averaged 
coronagraphic PSF's to the AO-only PSF. The most effective central image suppression is seen over 
a few diffraction rings outside the image plane stop for the intermediate value of J^. 
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Fig. 4. — The effect of increasing AO correction is seen here. The same 100 independent realizations 
of a Kolmogorov-spectrum phase screen representing atmospheric wavefront aberration corrected by 
AO systems with 18, 26, 36, 51 and 71 actuators across aD = 8.1m primary mirror, with D/vo = 30. 
The azimuthally averaged PSF's are normahzed to unity at their center. The Strehl ratios are 0.30, 
0.53, 0.69, 0.82 and 0.90 respectively. Note the widening plateaus of partial correction extending 
to the shoulder (solid triangle) at a distance of 9ao from the center of the image. Outside this 
plateau AO correction does not improve the image. 



-23- 




Fig. 5. — The effect of varying degrees of AO correction on increased dynamic range produced by a 
coronagraph with a AX/D occulting stop diameter. The dotted hues are azimuthally averaged, co- 
added PSF's for the same 100 reahzations of Kolmogorov-spectrum atmospheric phase disturbances 
with D/ro = 30. Sohd Hnes are renormahzed azimuthally averaged co-added PSF's of the same 
atmospheric phase realizations using a coronagraph with a Lyot stop fine-tuning parameter of 
J- = 0.5. The vertical scales arc identical, and the non-coronagraphic PSF's are normalized to 
unity at the origin. Renormalization is performed by dividing coronagraphic PSF's by the fractional 
throughput of the corresponding Lyot stops (see text). Panels a through f show the emergence of 
an annulus of image suppression outside the occulting stop edge as the Strehl ratio increases from 
53% to 90%. 
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Fig. 6. — Azimuthally averaged image profiles of observed PHARO _ftr-band data (solid) and and 
our simulations with D/ro = 10, and a parabolic AO high pass filter (dashed). Both images were 
normalized to contain unity power out to 4'.'0. 
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Fig. 7. — Azimuthally averaged PSF^s for the 1.2m Calypso telescope, with 7 actuators across the 
pupil diameter (~ 38 actuators within the pupil) and an occulting image stop of A\/D. (0'.'48 at 
A = 0.7;um). The Lyot stop fine tuning parameter T is 0.25 on the left, and 0.5 on the right (see 
text). The solid line shows the ^0-corrected PSF, the dashed line shows the image profile after 
the coronagraph. 
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Fig. 8. — Azimuthally averaged PSF^s for the 5m Palomar Hale telescope, with 16 actuators across 
the pupil diameter (~ 216 actuators within the pupil) and an occulting image stop of A\/D ((y/27 
at A = 2.2;um). The Lyot stop fine tuning parameter is 0.25 on the left, and 0.5 on the right (see 
text). The solid line shows the ^0-corrected PSF, the dashed line shows the image profile after 
the coronagraph. 
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Fig. 9. — Azimuthally averaged PSF^s for the 8.1m Gemini telescope, with 51 actuators across the 
pupil diameter (~ 2042 actuators within the pupil) and an occulting image stop of A\/D ((y.'16 at 
A = 1.6/um). The Lyot stop fine tuning parameter T is 0.25 on the left, and 0.5 on the right (see 
text). The solid line shows the 740-corrected PSF, the dashed line shows the image profile after 
the coronagraph. 
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Fig. 10. — Realizations of a star with a faint companion at 4A/D separation with D/vq = 15 for 
the Gemini 8.1m telescope assuming 2042 channels of AO. In (a) we show the primary star and a 
companion of Am = +5 with no coronagraph in place. Figures (b), (c) and (d) exhibit companions 
of Am = +5, +7.5 and +10 magnitudes respectively, with the primary occulted by an image plane 
stop of diameter AX/D (0^'16 in the iJ-band). 



